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Abstract 



The problem of kink stability of isothermal spherical self-similar flow in 
newtonian gravity is revisited. Using distribution theory we first develop a 
general formula of perturbations, linear or non-linear, which consists of three 
sets of differential equations, one in each side of the sonic line and the other 
along it. By solving the equations along the sonic line we find explicitly the 
spectrum, k, of the perturbations, whereby we obtain the stability criterion for 
the self-similar solutions. When the solutions are smoothly across the sonic 
line, our results reduce to those of Ori and Piran. To show such obtained 
perturbations can be matched to the ones in the regions outside the sonic 
line, we study the linear perturbations in the external region of the sonic line 
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(the ones in the internal region are identically zero), by taking the solutions 
obtained along the line as the boundary conditions. After properly imposing 
other boundary conditions at spatial infinity, we are able to show that linear 
perturbations, satisfying all the boundary conditions, exist and do not im- 
pose any additional conditions on L As a result, the complete treatment of 
perturbations in the whole spacetime does not alter the spectrum obtained 
by considering the perturbations only along the sonic line. 

PACS Numbers: 97.10.Bt, 97.60.-s, 98.35.Df, 98.62.Dm 

I. INTRODUCTION 

In this paper we re-consider the problem of kink instability of isothermal spherical self- 
similar flow in newtonian gravity. This problem was first considered by Ori and Piran [1] 
and found that a large class of solutions was not stable against kink perturbations. 

The kink modes result from the existence of sonic lines (points), along which disconti- 
nuities of (higher order) derivatives of some physical quantities can propagate. Such weak 
discontinuities cannot occur spontaneously. They are always the results of some singularities 
of the initial and/or boundary conditions. For example, they may occur due to the presence 
of angles on the surface of a body past which the flow takes place. In this case the first 
spatial derivatives of the velocity are discontinuous. They may also be formed when the cur- 
vature of the surface of the body is discontinuous (but without angles), in which the second 
spatial derivatives of the velocity are discontinuous. In addition, any singularity in the time 
variation of the flow also results in a non-steady weak discontinuity. The instability of such 
discontinuities is characterized by the divergence of them, and the blow-up may imply the 
formation of shock waves [2]. 

An example that discontinuities of derivatives can propagate along a sonic line is given 
by the linear perturbation, 5<f(r, x) = (pi(x)e kT , of a massless scalar field in 2 + 1 gravity, 
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which satisfies the equation [3], 

x(l-x) vV' + \ [(1 + 2k) - x(3 + 2k)] ^ - l -k Vl = /(*), (1.1) 

where a prime denotes the ordinary differentiation with respect to the indicated argument, 
f(x) is a smooth function of x, and x — 1 is the location of the sonic line 1 . From the above 
equation we can see that it is possible for ipi to have discontinuous derivatives across the 
sonic line x — 1. In fact, assume that <pi is continuous across x = 1 (but not its first-order 
derivative), we can write it in the form 

^i(x) = <pt(x)H{x - 1) + ^(x) [1 - H(x - 1)] , (1.2) 
where H(x — 1) denotes the Heavside (step) function, defined as 

ri, x>i, 

#(x-l) = (1.3) 
[0, x<\. 

Then, we find that 

tpx' = <pfH(x - 1) + <pf [1 - # (x - 1)] , 

^i" = ^"^(x - 1) + v?r" [i - # (* - 1)] + [yvr <Hx - 1), (i.4) 

where 8(x — 1) [= dH(x — l)/dx] denotes the Dirac delta function, and 

= Hm f^lMV (1.5) 

Lr J x^i+o \ dx J \ dx J 

Substituting Eq.(1.4) into Eq.(l.l) and considering the facts 

tf m (x) = tf(x), [1 - H(x)] m = [1 - tf(x)] , 

[1 - H(x)\ H(x) = 0, x5(x) = 0, (1.6) 



1 Note the different notations used here and those in [3]. In particular, here we use x in the place 
of z of [3]. 
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where m is an integer, one can see that Eq.(l.l) holds also on the sonic line x — 1 even 
when [yV]~ 7^ 0. This is because (1 — x) [</V]~ S(x — 1) = for any finite [</V]~- Thus, the 
quantity 



S( P, 



;i.7) 



has support only along the sonic line, where ( ) jX = d( )/dx. The above expression shows 
clearly how the discontinuity, 6(p x , of the first derivative of the perturbation Sip(r, x) 
propagates along the sonic line. When Re(k) > the perturbation grows exponentially as 
r — > 00, and is said unstable with respect to the kink perturbation. When Re(k) = 0, the 
perturbation oscillates in time r, and is usually also considered as unstable. When Re(k) < 
the perturbation decays exponentially and is said stable. Note that if the discontinuity 
happened on other places, say, 1 = i / 1, clearly Eq.(l.l) would not hold on x = x , 
because now (1 — x) [</V]~ S(x — xq) 7^ 0. This explains why the discontinuities are allowed 
only along sonic lines. 

From the above example we can also see that in general one may divide the whole space 
into three different regions, x > x c , x = x c , and x < x c , where the location of the 

sonic line. In each of the three regions, the perturbations, collectively denoted by Sy, satisfy 
a set of differential equations, 



G + (5y^\r,x,k) =0, (x>x c ), 
G~ (5y {m \r,x,k) =0, (x < x c ) , 



(1.8) 
(1.9) 

;i.io) 



where m — 0, 1, 2, and 



Sy(r, x) = Sy + (r, x)H (x - x c ) + 5y (r, x) [1 - H(x - x c )] , 

^> (Til)s ?v^), 



= lirn^ 



d m (5y + (r, x)) ± ^ d m (Sy-(r, x)) 



dx T 



dx T 



'1.11) 
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Therefore, in the presence of weak discontinuities the field equations actually have three 
different parts. Eq.(l.lO) describes the evolution of the discontinuities, Sy^ , along the 
sonic line, while Eqs.(1.8) and (1.9) describe the evolution of 8y ± in the regions x > x c and 
x < x c . Although 5y ± and 5y^ satisfy different equations, they have to match each 
other according to Eq.(l.ll). In addition, to have the perturbations physically acceptable, 
5y ± usually needs to satisfy some boundary conditions in the regions x > x c and x < x c , 
such as, the regularity conditions on the symmetry axis, and those at spacelike infinities. 

In [1] the field equation (1.10) along the sonic line was studied, by assuming that 
5y~(r,x) = in the region \x\ < \x c \. Hanawa and Matsumoto applied such a study to 
a fluid with a polytrope equation of state, p = Kp 1 , where p and p denote, respectively, the 
pressure and mass density of the fluid, and K and 7 are two constants [4]. Lately, Harada, 
and Harada and Maeda generalized such a study to the relativistic case [5,6]. Since Sy~(r, x) 
vanishes identically, it is necessary that 5y + (r,x) does not, in order to have non-vanishing 
§y(™) Then, a natural question is how the perturbations obtained along the sonic line 
match to the ones in \x\ > \x c \? 

To answer the above question, in this paper we first solve the equations along the sonic 
line for 5y c {r) = 5y(r,x c ), and find explicitly the spectrum, k, of the perturbations. In [1] 
only the case where the self-similar solutions are smooth across the sonic line was studied. 
Here we generalize such a study to include the case where the self-similar solutions can 
have discontinuous derivatives. Once 5y c {r) is found, we then solve the equations for the 
perturbations Sy + (r, x), by taking Sy c (r) and its first derivatives as the boundary conditions 
[cf. Eqs.(5.6) and (5.7)], 



Sy + (r,x) _ =5y c (r), (1.12) 

X — Xc 

6y+(T,x) _ =u(t). (1.13) 

X — Xc 



In addition to these, we also impose conditions at spatial infinity x = 00 [cf. Eq.(5.8)], 

Perturbations grow slower than the background solutions. (1-14) 
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In general, conditions (1.12) are already sufficient to determine all the integration constants. 
Thus, conditions (1.13) and (1.14) represent additional restrictions. By carefully analyzing 
the perturbations, we are able to show that there exist solutions to the perturbation equations 
that satisfy all the boundary conditions (1.12) - (1.14), an d do not alter the spectrum k, 
obtained by considering the perturbations only along the sonic line. 

Specifically, the paper is organized as follows: in Sec. II we review some main properties 
of the self-similar solutions, and pay particular attention on their asymptotic behavior near 
the sonic line and at the spatial infinity. In particular, we work out the equations that the 
discontinuities of the first-order derivatives of the self-similar solutions have to satisfy across 
the sonic line. In Sec. Ill, by using distribution theory we first develop a general formula 
of perturbations, which are valid for both linear and non-linear perturbations. It consists 
of three different sets, one in each side of the sonic line, and the other along it. Then, in 
Sec. IV we restrict ourselves to the perturbations along the sonic line and find explicitly 
the spectrum, k, of the perturbations. When the solutions are analytical across the line, we 
re-discover the results of Ori and Piran. When they are not, we find explicitly the stability 
criterion for those solutions in terms of x c . In Sec. V, we consider the linear perturbations 
Sy + (r, x) in the region \x\ > \x c \, by taking the solutions 5y c (r) and its derivatives obtained 
along the sonic line in Sec. IV as the boundary conditions. The paper is ended with Sec. 
VI, where our main conclusions are presented. 

II. FIELD EQUATIONS AND SELF- SIMILAR SOLUTIONS 

The spherical flow of self-gravitating isothermal perfect fluid is described by the hydro- 
dynamic equations, 



v, t + vv, r + 



cj,p tr Gm 



= 0, 



(2.1) 





(2.2) 



m t = — 47rr 2 pv, 



(2.3) 
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m r = 4:7ir 2 p, (2.4) 

where v >t = dv/dt, and r, t, c s , G, v(t,r), p(t,r) and m(t,r) denote, respectively, the radial 
coordinate, the time, the speed of sound, the gravitational constant, the velocity, the mass 
density, and the total mass enclosed within the radius r of the perfect fluid at the moment 
t, with the equation of state 

P = c 2 sP , (2.5) 

where p denotes the pressure of the fluid. The boundary conditions at r = are 

u(t,0) = 0, p(t,0) = 0. (2.6) 

Note that Eqs.(2.1) - (2.4) are unchanged under the coordinate transformations, 

t = i! + e , r = r', (2.7) 

where e is a constant. We call such a gauge as the Galilean gauge, and shall be back to it 
when we consider the linear perturbations in Sec. V. 

Following Ori and Piran, we first introduce the dimensionless quantities, 

x = — t = — In ( — ^ 

c s f \ t J ' 

V = -, D = —fr\ M=—, 2.8 

c s cj cjr 

where t is a dimensional constant with the dimension of length, which, without loss of 
generality, shall be set to one. Then, it can be shown that Eqs.(2.1)-(2.4) take the form 

(V -x)V x + E jX -V T + i(M - 2) = 0, (2.9) 

(V-x)E, x + V x -E :T = 0, (2.10) 

xM jX + M - e E = 0, (2.11) 

xM >x + M >T - -Ve E = 0, (2.12) 

Ob 

where E = ha(D). Since in this paper we are mainly concerned with the gravitational 
collapse of the fluid, from now on we shall restrict ourselves only to the region where t < 0, 



in which a singularity is usually developed at the moment t — 0. From Eq.(2.8) we can see 
that this also means x < 0. 

Self-similar solutions are given by 

F(t,x) = F m (x), (2.13) 

where F = {E, M, V}. Setting all the terms that are the derivatives of r be zero in 
Eqs.(2.9)-(2.12), we find that 

(V ss - x) VJ + EJ + - (M ss - 2) = 0, (2.14) 

x 

{V ss - x) EJ + V ss ' = 0, (2.15) 

xMj + M ss - e E °° = 0, (2.16) 

xM s J - -e E °°V S s = 0. (2.17) 
x 

The boundary conditions of Eq.(2.6) become 

V ss (0) = = D ss (0). (2.18) 

From Eqs.(2.16) and (2.17) we find that 

M ss = e E -(l-^y (2.19) 

The solutions of these equations have been studied extensively [7-10]. In particular, it 
was shown that they must pass a sonic line at x — x c , defined by 2 

V„(x c )-x c = l, (2.20) 

which divides the (t,r)-plane into two regions: the interior, \x\ < \x c \, where V 88 — x — 1 < 0, 
and the exterior, \x\ > \x c \, where V ss — x — 1 > 0, as shown in Fig. 1. The motion of 
inwardingly moving perturbations is given by dr/dt = v(t,r) — c s , or 



2 The line V ss (x c ) — x c = — 1 is also a sonic line, with the segment t < being saddle. However, 
this part is physically irrelevant, so in this paper we shall not consider it. 
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dx 
~dt 



1 



(V 9a -X-1) 



'<o, 


\x\ 


> 


\%c\ 


< =0, 


\x\ 




\%c\ 




\x\ 


< 





(2.21) 



where t < 0. Thus, in the exterior all the perturbations, including those that move in- 
ward with respect to the fluid, are dragged outward in terms of x, while in the interior 
perturbations can move in both directions. As a results, perturbations cannot penetrate 
from the exterior into interior [1]. this observation is important when we consider the kink 
perturbations. 
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FIG. 1. The (t, r)-plane. The line x = x c < represents the sonic line. The whole plane is 
divided into three regions: S7 + = {t, r : \x\ > \x c \}; £ = {t,r : x = x c }; and fl~ = {t, r : \x\ < \x c \}. 



Assuming that F ss (x) is at least C 2 in Q ± , it can be shown that the solutions near the 
sonic line x = x c are given by [8,10] 

V ss (x) = (1 + X c ) + 1 + — (X - X c ) - ( X - X cf + R V OG %c) , 

E ss (x) = In (-2x c ) - +Xc (x - x c ) + 1 (x - x c ) 2 + (x, x c ) , 

Ob q Ob ' £ 

M ss (x) = 2 - 2(1+Xc) (x - x c ) + (1 + Xc) 2 (2 + * c) (x - x c ) 2 + R M (x, x c ) , (2.22) 

x c x c 

for type 1, and 

V aa (X) ={1+X C )-—(X- X c ) - 2x 2( 2x + g) ( X ~ X c) + R E {X, X c ) , 

\ x 2 _ 2 

£ ss (x) = In (-2x c ) + — (x - x c ) + 2x 2(2x + 3) ^ x ~ Xc ^ + Re x< ^ ' 

M ss (x) = 2- 2(1 + Xc) (x - x c ) + 2 ^ (x - Xcf + Rm (x, x c ) , (2.23) 

X c X c 

for type 2, where -Rf #c) denotes the errors of the corresponding expansion. Note that 
when x c = —2 the two types coincide. For our present purpose, we can consider this 
degenerate case as a particular one of any of the two types. 

The qualitative behavior of the solutions near the critical point depends on the values 
of x c : (a) When |x c | < 1, the critical point is a saddle, and only two solutions pass through 
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it, one for each type, (b) When \x c \ > 1, the critical point is a node [9]. In this case, the 
directions that the solutions pass the critical point are further classified into primary and 
secondary directions. In particular, one of the directions, type 1 solutions for 1 < \x c \ < 2 
and type 2 solutions for \x c \ > 2, is called the secondary direction, while the other the 
primary direction. Only one solution can pass a node in the secondary direction, but an 
infinite number of solutions pass it along the primary direction [10]. Physically relevant 
solutions are these across the sonic line through nodes [cf. Fig. 2]. Thus, in the following 
we shall consider only the case where x c < —1. 




X 

FIG. 2. The (V ss , x)-plane. The region bound by the lines V ss (x) = x and x = is excluded, 
because in this region the mass of the fluid is negative, m < 0. The lines with arrows represent a 
possible trajectory of the fluid. 

It is possible to match solutions with different types across the critical point. But in this 
case the matching is only C°, that is, F ss (x) is only continuous, and its first-order derivatives 
normal to the surface x = x c are usually not [9]. Ori and Piran considered only the case 
where the matching is smooth, that is, solutions that belong to the same type in both sides 
of the sonic line. In this paper, we shall consider all the possibilities. 

Dividing the whole plane into three regions, and S, where Q + = {t,r : \x\ > \x c \}, 
S = {t,r : x — x c }, and Q~ = {t,r : \x\ < \x c \}, we can write F ss (x) as, 
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F S s(x) 



(F+(x), xeQ+, 

F c , x e s, 
F~(x), x e Q~. 



[2.24) 



Since F ss is C° across the surface x = x c , we must have 

E+(x c ) =E-(x c ) = ln(-2x c ), 
M+ (x c ) = M~ (x c ) = 2, 
V+ (x c ) = V~ (x c ) = l + x c , 



(2.25) 



as we can see from Eqs.(2.22) and (2.23), this is indeed the case. On the other hand, from 
Eqs. (2.19), (2.22) and (2.23) we find that M ss ' is also continuous and given by 

Mt: (*c) = M~' (x c ) = - 2 J±±^A. (2.26) 

x c 

Taking the limits x — > xf° in Eqs. (2. 14) - Eqs. (2. 17) and then subtracting them, we find 
that 

[VssT + [ E ssT = °> (x = x c ). (2.27) 

Note that in writing Eq.(2.27) we had used Eqs. (2.20) and (2.25). 

As assumed previously, F ss (x) is at least C 2 in fl ± , we can take derivative of Eqs. (2. 14) 
- (2.17) with respect to x, and obtain 



(V ss - x) V ss " + E ss " + (VJ - 1) VJ + -MJ - -i (M ss - 2) = 0, 



x 



(V ss - x) E ss " + V ss " + {Vss' - 1) Es/ = 0, 
xM ss " + 2M SS ' - e E °°Ess' = 0, 

xMs/' + Ms/ - -e E - (VJ + VsEss') + \e E °°V ss = 0, 



X 



(2.28) 
(2.29) 
(2.30) 
(2.31) 



where the quantities F ss (x) in the above equations should be understood as F+(x) in Q + 
and F~(x) in tt~ . The subtraction of the limits x — > xf° of Eqs.(2.28) - Eqs.(2.31) yields 

[Vss'T + [EssT - [VssT (l - [VssT) = 0, (2.32) 

[VssT + [EssT - [EssT + \ (\E s lt [VssT + [EssT [VssT) = o, (2.33) 



[M ss "] +2 [EsJ] =0, (x = x c ), 



(2.34) 
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where 



„+ ,. idF+(x)\ (dF-(x) 



[FJV = lim 'f w + lim [ ss 7 v 7 . (2.35) 

Note that in writing the above equations we had used Eqs.(2.25) and (2.27), so that 
Eqs.(2.30) and (2.31) yield the same equation, (2.34). From Eqs.(2.32) and (2.33), on the 
other hand, we find that 

[V SS T ([EJ} + - 3 [VJ] + + 4) = 0, (x = x c ), (2.36) 

which has the solutions, 

(a) [V sa T = 0, (2.37) 

(b) [Ej] + -3[V ss t + 4 = 0. (2.38) 

Clearly, Case (a) requires that the solutions inside and outside the surface x = x c belong to 
the same type, and in this case from Eqs.(2.22), (2.23) and (2.27) we find that 

[^ / r-3[K;r+4=i^Mx( 1, Type1, (^-=0). (2. 39) 

Fcl [-1, Type 2, V ' 

III. GENERAL FORMULA OF PERTURBATIONS 

In this section using distribution theory, we shall develop formulas that are valid for any 
kind of perturbations, linear or non-linear, and write down the equations of perturbations 
in and E separately. 

For any given function F(x, r) that is C 2 in ^ and C° across the surface be 
written as 3 , 

F(x, r) = F + (x, t)[1-H(x- x c )\ + F~(x, r)H (x - x c ) , (3.1) 
where F + (F~) is defined in Q + and 



3 Note that in Jl + (Q ) we now have x < x c (x > x c ). 
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\im o F + (x,r) =Km o F-(x,T) = F c (r). (3.2) 

In the present case, it is clear that 

F ± (x,t)=F±(x)+5F ± (x,t), (3.3) 

where F^(x) denotes the self-similar solutions of Eqs.(2.14) - (2.17) with the boundary con- 
ditions of Eqs. (2.18), (2.25) and (2.26). The perturbations SF ± (x, r) satisfy the conditions, 

\im + JF + (x,T) = lim q 6F-(x,t) = 8F c (t). (3.4) 

X ^X q X ^X q 

Then, using Eqs. (1.6), (3.1), (3.2), and the fact that 

H(x- x c ) 5 (x - x c ) = ^5 (x - x c ) , [1 - H (x - x c )] S(x- x c ) = ^8(x- x c ) , (3.5) 
we find 

F x = F+, x [1-H(x- x c )\ + F~ jX H (x - x c ) , 

F iXX = F + :XX [1 - H (x - x c )] + F~ iXX H (x - x c ) - [F :X ]~ 5 (x - x c ) , 

F (n) = F +W ^ _ H ^ _ + F -(n) H ^ _ ^ _ ^ [ F (n-i-l)j " ^ _ ^ (3 g) 



i=0 



where 





= lim 




X *x c \ 



d n F+ (x,t) \ ( d n F~ (x,t) \ 

dx" ) ~ V /' 

x(n) I \ — { x ~ x c) 

5 {X ~ Xc) = dx" ' 

/oo 
F(x, r)5 ik) (x - x c ) dx = (-l) k F^ (x c , r) . (3.7) 
-oo 

Inserting Eqs.(3.1)-(3.6) into Eqs. (2.9) - (2.12), and considering Eq.(3.5), we find that 
(V ss - x) 5V, X + V„ tX 8V + 5V5V, X - 5V, r + 5E, X + -5M = 0, (3.8) 

Ob 

(V ss - x) 5E, X + E SS JV + 5V5E, X - 5E, T + 5V >X = 0, (3.9) 

x5M jX + 5M - e Eas (e 5E - l) = 0, (3.10) 

x5M tX + 5M )T + -e Ess [v ss (l - e 5E ) - e SE 8V\ = 0. (3.11) 

Ob 
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Again, the quantities F, 5F in the above equations should be understood as 
F + , 5F+ (F~, 6F~) in fi+ The boundary conditions (2.6), together with those 

of Eq.(2.18) require 



6V-(t,0) = = SD-(t,0). 

Subtracting the limits x xf° of Eqs.(3.8) - (3.11), on the other hand, we find, 

(1 + 5V C ) W c (t) + U c (r) + [VsJf W c = 0, 
(1 + 5V C ) U c (t) + W c (r) + [Ejf 5V C = 0, 
[SM, x r = 0, (x = x c ), 



(3.12) 



(3.13) 
(3.14) 
(3.15) 



where 



U c (r) = [5E /x ]- , W c (r) = [SV,. 



(3.16) 



and [SF >X \ is defined as 



± lim 



'9(5F- (x,r))' 



9a; / z^xj V dx 
Note that in writing Eqs.(3.13) -(3.15), we used the fact that 



(3.17) 



|Uf+(.,t) 



= £-5F~ (x,r) 
„+o ar 



(3.18) 



In addition, Eqs.(3.10) and (3.11) give the same equation, (3.15). Combining Eq.(2.27) with 
Eqs.(3.13)-(3.15) we find the following solutions, 



(i) 5V C = 0, U c (t) = -W c (t), [5M iX ]' = 0, 

(ii) 6V c = -2, U c (t) = W c (r) + 2 [V sa 'r , [SM tX ]~ = 0, 

(iii) U c (t) = -W c (t) = [VJ]- , [5M,J^ = 0, (x = x c ). 



(3.19) 
(3.20) 
(3.21) 



On the other hand, taking the derivatives of Eqs.(2.9) - (2.12) with respect to x, we find 



that 
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(V - x) V xx + E, xx - V XT + (V x -1)V X + -M, x — \ (M — 2) = 0, (3.22) 



x ' x 2 



{V - x) E iXX + V xx - E, XT + (V x - 1) E iX = 0, (3.23) 

xM )XX + 2M X - e E E, x = 0, (3.24) 

xM, xx + M iTX + M x - -e E (V x + VE, X ) + \ve E = 0. (3.25) 

x x 2 

Substituting Eqs.(3.3) - (3.6) into the these equations and using Eq.(3.5), we find that each 
of these equations can be written in the form, 

G + (t, x)[l-H{x- x c )\ + G~(t, x)H (x - x c ) + G c (r) 5 (x - x c ) = 0, (3.26) 

which is equivalent to 

G + (t,x)=0, (xGfi+) (3.27) 

G-(t,x)=0, (iGfij (3.28) 

G c (r) = 0, (igE). (3.29) 

It can be shown that Eqs.(3.27) and (3.28) take the same forms as these given by Eqs.(3.22) 
- (3.25) after replying F by F ± . But, on the surface x = x c Eq.(3.29) gives 

[V,x]~ + [E,x]~ — 0, (3.30) 
[M x ]- = 0. (3.31) 

Eq.(3.31) is consistent with Eqs.(2.26) and (3.15), while Eq.(3.30) together with Eq.(2.27) 
yields 

U c (t) + W c (t) = 0. (3.32) 

On the other hand, subtracting the limits x — > xf° of Eqs.(3.22)-(3.25), we find that 

_ (i + SVc) [Sv ^y _ [SEtXX }- - 5V C [VJT 
+ (l - [V ss f - [5V, X } + ) W c (r) - [V ss '\- [5V, X } + = 0, (3.33) 
<EM _ (i + SVc) [5 E tXX }- - [5V, XX \~ - 5V C [EJT 
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+\ (2 - [V ss '] + - [5VJ + ) U c (r) - \ ([Ej} + + [6EJ + ) W e (r) 

- l - ([Vjf [5E, X ] + + [Ejf [5V, X ] + ) = 0, (3.34) 

x c [5M >XX ]- - 2x c (l - e 5E <) [Ejf + 2x c e &Ec U c {r) = 0, (3.35) 

([Ejf + U c (t)) 5V C = 0. (3.36) 

Equations (3.8) - (3.11), (3.19) - (3.21), and (3.32) - (3.36) consist of the full set of differential 
equations for the perturbations 5F ± {t, x) in the whole spacetime, which consists of three 
regions, Vt^ and E. Note that these equations are exact, and so far no approximation has 
been made. In Sec. V we shall consider their linearized perturbations. 



IV. KINK STABILITY 

Kink stability is the study of the perturbations of Eqs.(3.19) - (3.21), and (3.32) - (3.36) 
along the critical line x = x c . To solve these equations for 5F ± (t, x), in addition to (3.12), 
other boundary conditions need to be given. Following Ori and Piran, we shall impose the 
following conditions: Assume that the perturbations turn on at the moment t = t or r = 
[cf. Fig. 1], then we require 

(A) the perturbations initially vanish in the interior, 

5F~(t = 0,x) = 0, x E tt~, (4.1) 

(B) the perturbations be continuous everywhere, and in particular across the sonic line, 

[<SF]- = 0, (x = x c ), (4.2) 

(C) 8E ± :X and 5V ± jX be discontinuous at the sonic line, 

U c (t) = [SE >x ]-^0, W c (r) = [5V X ]- ^ 0, (x = x c ) . (4.3) 

From the above we first note that Eq.(4.1) remains true for all the r > 0, because the 
perturbations cannot penetrate the sonic line from exterior into the interior, as we noticed 
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in Sec. II, and 5F~(x,t) = are indeed solutions of Eqs.(3.8) - (3.11) in f2~. Thus, in the 
following we need to consider only the perturbations in Q + as well as those along the sonic 
line x = x c . 

The studies of the perturbations in Q + will be considered in the next section, while in 
the rest of this section we shall consider only Eqs.(3.19) - (3.21) and (3.32) - (3.36). As first 
noticed by Ori and Piran, the evolutions of these two sets of equations are separable. In 
particular, Conditions (4.1) - (4.3) together with Eqs.(3.19) - (3.21) and (3.32) - (3.36) are 
already sufficient to determine the evolution of [<5-F\rp uniquely. To show this, let us first 
notice that Conditions (A) and (B) imply 



5F c (t) =8F- (t,x c )=0, 
[oF :X \ = ■ 



[5F, X } + - (4-4) 



dx 

Then, we can see that only the case of Eq.(3.19) is consistent with these expressions. Thus, 
in the following we shall discard the cases described by Eqs.(3.20) and (3.21). Hence, one 
can see that Eq.(3.36) is satisfied identically, while from Eqs.(3.33) - (3.35) we find 



dWc 
dr 



+ - (4 + [EJ} + - 3 [VJ] + - 4 [VJ]-) W c - W 2 C = 0, (4.5) 



where in writing this equation we had used Eq.(3.32). To study it, let us consider the two 
cases defined by Eqs.(2.37) - (2.38) separately. 

Case A) [V ss '] = 0: In this case the internal and external solutions belong to the same 
type, and we have 



F£(x e ) = F- S '(x c ) = FJ, (4.6) 



for which Eq.(4.5) becomes 



^£ = kW c + Wl k = (2V c '-l). (4.7) 
dr 

This is exactly the equation obtained by Ori and Piran [1]. In particular, it was shown that 
when V c ' > 1/2 the self-similar solutions are unstable against the kink perturbations. 
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W c {r) ~ W Q e 
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Case B) \V SS '\ ^ 0: This is the case of Eq.(2.38), from which we find that Eq.(4.5) 
takes the same form as that of Eq.(4.7), but now with 

k=[VjV. (4.8) 

The dynamics of small perturbations are determined by the linear part of Eq.(4.7), 

^ kW c , (4.9) 
dr 

which has the solution, 

( ±oo, k > 0, 

W , k — 0, (as t -> oo), (4.10) 
10, fc < 0, 

where Wq is an integration constant and denotes the initial profile of the perturbation. 
Therefore, when k > the corresponding self-similar solutions are unstable against kink 
perturbations, and when k < 0, they are stable. To find out the sign of k we need to 
consider the two cases where the external solution is Type 1 or 2 separately. When it is type 
1, then the internal solution must be type 2, and then from Eqs.(2.22) and (2.23) we find 

1x1-2 f > 0» x c < -2, 
k = V+\x e ) - V-'(x c ) = l -^ T - = (4.11) 

Fc| i<0, -l>x c >-2. 

Thus, for this kind of matching, the self-similar solutions are stable for — 1 > x c > —2, 

and unstable for x c < — 2. When the external solution is type 2, the internal solution must 

be type 1. Then, the sign of k is just opposite to the last case, and we conclude that 

the corresponding self-similar solutions now are unstable for —1 > x c > —2 and stable for 

x c < -2. 

It is interesting to notice that Eq.(4.10) allows the following general solution [1], 

which can diverge for certain choice of the initial profile Wo- Ori and Piran argued that this 
divergence could be the indication of the formation of a shock wave. 
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V. LINEAR PERTURBATIONS IN REGION U+ 



In this section, we shall study the linear perturbations of the self-similar solutions given 
in Sec. Ill in the region Q + , where \x\ > \x c \, with SF + (t,x c ) and SF + ^ X (t,x c ) obtained in 
the last section as the boundary conditions. Let us first write the perturbations SF + (r, x) 

as 

5F + (r,x) = eF 1 (x)e k \ (5.1) 

where e is a very small quantity, and F\ (x) denotes perturbations. It is understood that there 
may have many perturbation modes for different k. Then, the general linear perturbation 
will be the sum of these individual ones. 

To the first order of e, the substitution of Eq.(5.1) into Eqs.(3.8) - (3.11) yields 

(V ss - x) V,' + Ei + (VJ -k)Vi + -Mi = 0, (5.2) 

Ob 

(V ss - x) Ei' + Vi + E ss 'Vi - kEi = 0, (5.3) 
xMi' + Mi — e Ea °Ei = 0, (5.4) 

xMi' + kMi - -e Esa (V ss Ei + Vx) = 0, ( x e Q + ) , (5.5) 

x v ' 

where, without causing any confusions, in writing the above equations we had dropped the 
superscript "+". 

Physically acceptable solutions of the above equations must satisfy some boundary con- 
ditions at the sonic line X — 0C c 

well as at the spatial infinity x = oo. At x — x c , the 
boundary conditions are the solutions obtained in the last section along x = x c . From 
Eqs.(3.15), (3.32), (4.4) and (4.10) we find that these conditions read 

Vi (x c ) = Ei (x c ) = Mi (x c ) = 0, (5.6) 
V 1 '(x c ) = -Ei'(x c ) = W , M 1 '(x c ) = 0, (x = x c ), (5.7) 

where Wo is the integration constant given in Eq.(4.10). In the neighborhood of x = x c , the 
functions F ss (x) behave as those given by Eqs.(2.22) and (2.23). 
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On the other hand, at the spatial infinity x = oo (or r = oo), we require that the 
perturbations Fi(x) should not grow faster than the background solutions F ss (x), that is, 

Fi(x) grows slower than F ss (x), as x — > oo, (5.8) 

where Fi = {Ei, Vi, Mi}. The asymptotic behavior of F ss (x) is given by [9], 

VM) = Vx - <^> + % + + 2,(^-6)1 + Q , 4) 

x x^ ox- 3 v y 

F ss (x) = In (Do,) - (jP ^7 2) + O (x- 4 ) , (5.9) 

as x — > oo, where Doo and Voo are two arbitrary constants. 

The asymptotic behavior of M ss (x) and x — > oo can be obtained from 

Eqs.(2.19), (2.22), (2.23) and (5.9). 

Since F ss (x) is well-defined in the region x e (— oo, — |x c |), one can see that there always 
exist solutions to Eqs.(5.2)- (5.4), which in general has three free parameters. However, 
Eq.(5.5) represents a constraint, and shall reduce the number of the free parameters from 
three to two. On the other hand, one can show that Mi(x c ) = holds, provided that 
Vi(x c ) = Ei(x c ) = is true. Thus, the boundary conditions given by Eq.(5.6) are already 
sufficient to determine the two free parameters uniquely. Then, an important question rises: 
Do Eqs.(5.7) and (5.8) give further restrictions on these parameters? In particular, once all 
the conditions of Eqs.(5.6) - (5.8) are taken into account, do the values of the spectrum k 
still remain the same as those given by Eq.(4.7) or (4.8)? If not, is it possible that Re(k) 
cannot be positive any more, and that all the unstable modes found in the last section are 
due to the incomplete treatment of the problem? 

Before addressing these questions, we first consider the problem concerning gauge modes. 
Under the coordinate transformations (2.7), we have 

t — > t + ee T , x — * x (1 + ee T ) , 

F ss (x) -> F ss (x) + tF[ 9 \x)e\ (5.10) 

where 
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Fi 9 \x) = xF' ss (x). 



(5.11) 



Comparing it with Eq.(5.1) we find that the gauge mode corresponds to the case where 
k — 1. This mode is not physical, and completely due to the gauge transformations (2.7). 
Thus, in the following we shall discard this case and consider only the case where k ^ 1. 
Then, from Eqs.(5.4) and (5.5) we find 

1 



Mi 



(1 - k)x 

Combining it with Eqs.(5.2) and (5.3) we find 

1 



[(V ss -x) E 1 + V 1 ]e 



(5.12) 



Eh! 



(AV X + BE,) , 



where 



(V ss -x) z -l 

Vi'= n/ \ 2 (CV 1 + DE 1 ), 
(V ss - x) - 1 



A = (Ves ~ k) - (V ss ~ X) E ss ' 



(5.13) 
(5.14) 



„E 3 



(l-k)x 2 



B = (Vss ~ x) 



k 



[1 - k)x 2 



C = E ss ' - (V ss - x) 

\2 



(V ss ' -k) + 



1 - k)x 2 



D = 



^- x K^-k. 



(5.15) 



(1 - k)x 2 

In the following we shall consider the asymptotic behavior of E\ and V\ at x — x c and x = oo 
separately. First, in the neighborhood of x c we expand them as 



E 1 = E% + El (x - x c ) + R l E (x, x c ), 
V l = V? + V, 1 (x - x c ) + R\,(x, x c ). 

Substituting the above expressions into Eqs.(5.13) and (5.14) we find that 

[l-k)[(k-2)x c -2}-2. 



(5.16) 



k(l-k)x c + 2 

-v; 1 + kEl + 



■v?, 



" x 1 1 l+Xc V*, (Typel), 



Xr 



(5.17) 
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for type 1 solutions, and 

= k[(l-k)x c -2] Q 

1 k(l-k)x c + 2 

E{ = -V, 1 + kE[ - -V?, (Type 2), (5.18) 

x c 

for type 2 solutions, where and V^ 1 are arbitrary constants in both types. Applying the 
boundary conditions (5.6) and (5.7) to the above solutions we find that 

£0 = y0 = , El — -Vj 1 = —W , (5.19) 

for a given k. 

On the other hand, setting e = 1/rr we can expand F 1 at the neighborhood of e = as, 

E' = e° 1 + e\e + 0(e 2 ), 
M 1 = m\ + m}e + O (e 2 ) , 

y 1 = vl + vie + O (e 2 ) . (5.20) 

Inserting these expressions and the ones given by (5.9) into Eqs.(5.2)-(5.5) we find that all 
the coefficients, /{, are zero, unless k — 0, for which we have 



Doo (e° - v°) , 
e} = 0, = -Dooe?, 

^oo (^oo - V^e?) e°, (Type 1), (5.21) 



for type 1, and 



^Doo (e? - ) , 
e\ = 0, v{ = -D oa e\, 

D oo (D oo -V oo e 1 )e° 1 , (Type 2), (5.22) 



for type 2, where e\ and i>J are arbitrary. However, in each of the cases conditions (5.8) are 
satisfied for any given k. 
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In review all the above, we conclude that the considerations of the perturbations outside 
the sonic line x = x c and their matching to the ones obtained along the line do not impose any 
additional conditions on the spectrum, k, of the kink perturbations, obtained by considering 
the linear perturbations only along the sonic line. 

VI. DISCUSSIONS AND CONCLUDING REMARKS 

In this paper we first gave a general review over the weak discontinuities of a fluid along 
a sonic line by providing a concrete example, and showed how such a discontinuity can 
propagate along the line. Then, in Sec. II we summarized the main properties of the self- 
similar solutions of an isothermal spherical self-similar flow. In particular, we investigated 
the discontinuities of these solutions across the sonic line, which is useful and necessary for 
the studies of the kink stability to be considered in the following sections. In Sec. Ill, using 
distribution theory we first developed a general formula of perturbations in the background 
of a self-similar solution. For such a development, we assumed that the background is at 
least C 2 outside of the sonic line and C° across it. The field equations were divided into 
three different sets, one in the region x > x c , one in the region x < x c , and the other along 
the sonic line x = x c . Any perturbations must satisfy all of these equations. In Sec. IV, 
following Ori and Piran [1] we considered the differential equations along the sonic line and 
found their explicit solutions. When the background is smooth across the sonic line, we 
re-derived the results of Ori and Piran. When the background is not smooth, we found the 
stability criterion for the self-similar solutions. In particular, if the external solution is type 
1 and the internal is type 2, the solution is stable for — 2 < x c < —1, and not stable for 
x c < —2. If the external solution is type 2 and the internal is type 1, it is just the other way 
around. That is, it is stable for x c < —2, and not stable for —2<x c <—l. 

Up to this point, it was not clear whether the solutions of the perturbations obtained 
along the sonic line can be matched to the ones outside the sonic line. Since the perturbations 
in the region |x| < \x c \ is assumed to be zero identically, the ones in the region |x| > \x c \ are 
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necessarily non-zero. Otherwise, it is impossible to get a non-zero discontinuity across the 
sonic line. To answer the above question, we studied the linear perturbations in the region 
\x\ > \x c \, by taking the solutions obtained along the sonic line as the boundary conditions. 
By properly imposing other boundary conditions at the spatial infinity, we were able to show 
that there always exist linear perturbations in the region \x\ > \x c \, which satisfy all the 
boundary conditions for a given k, where k is the spectrum of the perturbations obtained 
along the sonic line. As a result, the studies of perturbations along the sonic line carried out 
before are consistent with the complete treatment of the problem in the whole spacetime. 
In particular, the spectrum of perturbations obtained along the sonic line remains the same 
even after the perturbations in the whole spacetimes are considered. 

Finally, we note that it would be very interesting to generalize the above studies to the 
relativistic subject that is under our current investigations. 
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